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Halo phenomena in deformed nuclei are investigated within a deformed relativistic Hartree Bo- 
goliubov (DRHB) theory. These weakly bound quantum systems present interesting examples for 
the study of the interdependence between the deformation of the core and the particles in the halo. 
Contributions of the halo, deformation effects, and large spatial extensions of these systems are 
described in a fully self-consistent way by the DRHB equations in a spherical Woods-Saxon basis 
with the proper asymptotic behavior at large distance from the nuclear center. Magnesium and neon 
isotopes are studied and detailed results are presented for the deformed neutron-rich and weakly 
bound nucleus ““Mg. The core of this nucleus is prolate, but the halo has a slightly oblate shape. 


This indicates a decoupling of the halo orbitals from the deformation of the core. 


The generic 


conditions for the occurence of this decoupling effects are discussed. 


PACS numbers: 21.10.Gv, 21.10.-k, 21.60.Jz, 21.60.-n 


The “shape” provides an intuitive understanding of 
spatial density distributions in quantum many-body sys- 
tems, such as molecules [1], atoms [2], atomic nuclei [3], 
or mesons [4]. Quadrupole deformations play an im- 
portant role in this context. The interplay between 
quadrupole deformation and the weak binding can re- 
sult in new phenomena, such as “quadrupole-bound” an- 
ions [5]. 

Halo phenomena in nuclei are driving forces for the 
development of the physics with radioactive ion beams. 
They are threshold effects [6] and have been first ob- 
served in the weakly bound system !!Li [7]. Considering 
that most open shell nuclei are deformed, the interplay 
between deformation and weak binding raises interesting 
questions, such as whether or not there exist halos in 
deformed nuclei and, if yes, what are their new features. 

Calculations in a deformed single-particle model [8] 
have shown that valence particles in specific orbitals with 
low projection of the angular momentum on the symme- 
try axis, can give rise to halo structures in the limit of 
weak binding. The deformation of the halo is in this 
case solely determined by the intrinsic structure of the 
weakly bound orbitals. Indeed, halos in deformed nu- 
clei were investigated in several mean field calculations 
in the past [9-11]. However, in Ref. [12], it has been 
concluded that in the neutron orbitals of an axially de- 
formed Woods-Saxon potential the lowest- component 
becomes dominant at large distances from the origin and 
therefore all Q7 = 1/2* levels do not contribute to de- 
formation for binding energies close to zero. Such argu- 
ments raise doubt about the existence of deformed halos. 
In addition, a three-body model study [13] suggests that 
it is unlikely to find halos in deformed drip line nuclei 
because the correlations between the nucleons and those 


due to static or dynamic deformations of the core inhibit 
the formation of halos. 

Therefore a model which provides an adequate de- 
scription of halos in deformed nuclei must include in a 
self-consistent way the continuum, deformation effects, 
large spatial distributions, and the coupling among all 
these features. In addition it should be free of adjustable 
parameters that make predictions unreliable. Density 
functional theory fulfills all these requirements. Spher- 
ical nuclei with halos have been described in the past 
successfully in this way by the solution of either the 
non-relativistic Hartree-Fock-Bogoliubov (HFB) [14-16] 
or the relativistic Hartree Bogoliubov (RHB) equations 
[17-19] in coordinate (r) space. However, for deformed 
nuclei the solution of HFB or RHB equations in r space is 
a numerically very demanding task. In the past consider- 
able effort has been made to develop mean field models ei- 
ther in r space or in a basis with an improved asymptotic 
behavior at large distances |11, 20-25]. In particular, an 
expansion in a Woods-Saxon (WS) basis was shown to 
be fully equivalent to calculations in r space [23]. 

In the present investigation, we therefore study halo 
phenomena in deformed exotic nuclei within a DRHB 
model using a spherical WS basis. The RHB equations 
for the nucleons read [17, 26] 
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where Ep is the quasiparticle energy, A the chemical po- 
tential, and hp is the Dirac Hamiltonian [27-31] 


hp =a-p+V(r)+B(M + S(r)). (2) 


Neglecting here for simplicity spin and isospin degrees of 
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FIG. 1: (Color online) Density distributions of **Mg with the 
z-axis as symmetry axis: (a) the proton density (for x < 0) 
and the neutron density (for x > 0), (b) the density of the 
neutron core, and (c) the density of the neutron halo. In each 
plot, a dotted circle is drawn for guiding the eye. 


freedom the pairing potential reads, 
A(r1, r2) = V”? (r1, r2)K(r1, r2), (3) 


with a density dependent force of zero range in the 
particle-particle channel 


VP? (r1, r2) = Vo ô(rı — r2) (1 -= eeo) Lisp (4) 


and the pairing tensor « which is defined in the conven- 
tional way. 

For axially deformed nuclei with spatial reflection sym- 
metry, we represent the potentials and densities in terms 
of the Legendre polynomials, 


=X fr(r)Py(cos0), A=0,2,4,:. (5) 
à 


For fixed quantum numbers Q” the deformed quasipar- 
ticle wave functions U(r) and V(r) in Eq. (1) are ex- 
panded in a spherical WS basis (for details see Ref. [23]). 

The calculations are based on the density functional 
NL3 [32] and the pp interaction (4) with the parameters 
Psat = 0.152 fm~?, Vo = 380 MeV-fm?, and a cut-off 
energy E2? = 60 MeV in the quasi-particle space. These 
parameters reproduce the proton pairing energy of the 
spherical nucleus ??Mg obtained from a spherical RHB 
calculation with the Gogny force D1S. A spherical box 
of the size Rmax = 20 fm and the mesh size Ar = 0.1 
fm are used for generating the spherical Dirac WS basis 
of Ref. [23] which consists of states with j < 2th. An 
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FIG. 2: (Color online) Single neutron levels with the quan- 
tum numbers Q” around the chemical potential (dotted line) 
in the canonical basis for “4Mg as a function of the occupation 
probability v?. The order n, 27, and the main WS compo- 
nents for orbitals close to the threshold are also given. The 
dashed line corresponds to the BCS-formula with an average 
pairing gap. 


energy cutoff E$ = 100 MeV is applied to truncate the 
positive energy states i in the WS basis and the number of 
negative energy states in the Dirac sea is taken to be the 
same as that of positive energy states in each (£, j)-block. 


In the present study of Mg isotopes, the last nucleus 
within the neutron drip-line is **Mg. Of course, it is 
difficult to predict the position of the drip-line precisely 
for nuclei so far from the experimentally known area and 
therefore the results discussed in the following have to be 
taken as generic results. In this study *°Mg is an almost 
spherical nucleus. The neighboring nucleus “4M¢g is well 
deformed (62 = 0.32) and weakly bound with the two- 
neutron separation energy Sən = 0.44 MeV. Therefore 
this nucleus is taken here as an example for a detailed 
investigation. The density distributions of all protons 
and all neutrons in “4Mg are shown in Fig. la. Due to 
the large neutron excess, the neutron density not only 
extends much farther in space but it also shows a halo 
structure. The neutron density is decomposed into the 
contribution of the core in Fig. 1b and that of the halo 
in Fig. lc. Details of this decomposition are given fur- 
ther down. We find that the core of 44Mg is prolate, and 
that the halo has a slightly oblate deformation. This in- 
dicates the decoupling between the deformations of core 
and halo. 

Weakly bound orbitals or those embedded in the con- 
tinuum play a crucial role in the formation of a nuclear 
halo. For an intuitive understanding of the single par- 
ticle structure we keep in mind that HB-functions can 
be represented by BCS-functions in the canonical basis 
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and show in Fig. 2 the corresponding single neutron spec- 
trum. As discussed in Ref. [33] the single particle energies 
in the canonical basis €p = (k|hp|k) shown in Fig. 2 are 
expectation values of the Dirac Hamiltonian (2) for the 
eigenstates |k} of the single particle density matrix ô with 
the eigenvalues vz. The spectrum of ô has a discrete part 
with v? > 0 and a continuous part with v? = 0. Obvi- 
ously only the first part contributes to the HB-wave func- 
tion and only this part is plotted in Fig. 2. This part of 
the spectrum €x is discrete even for the levels in the con- 
tinuum. Of course, this is only possible because the wave 
functions |k} are not eigenfunctions of the Hamiltonian. 
As long as the chemical potential An is negative, the cor- 
responding density p(r) is localized [15] and the particles 
occupying the levels in the continuum are bound. 

The orbitals in Fig. 2 are labeled by the conserved 
quantum numbers Q and m. The character n numbers 
the different orbitals appearing in this figure accord- 
ing to their energies. The neutron Fermi energy lies 
within the pf shell and most of the single particle lev- 
els have negative parities. Since the chemical potential 
An = —230 keV is relatively small, orbitals above the 
threshold have noticeable occupation due to pairing cor- 
relations. For example, the occupation probabilities of 
the 5th (Q7 = 7/27) and the 6th (Q7 = 1/27) orbitals 
are 27.2% and 14.3%. 

As we see in Fig. 2 there is a considerable gap between 
the two levels with the numbers n = 2 and n = 3. The 
levels with €can < —2.5 MeV contribute to the “core”, 
and the other remaining weakly bound and continuum 
orbitals with €c¢an > —1 MeV naturally form the “halo”. 
Therefore we decompose the neutron density p"(r) into 
two parts, one part coming from the orbitals with canon- 
ical single particle energies Ecan < —2.5 MeV (called 
“core” ) and the other from the remaining weakly bound 
and continuum orbitals (called “halo”). The spherical 
components of these densities [i.e. the contribution of 
A = 0 in Eq. (5)] are plotted together with that of the 
total neutron density in Fig. 3a. It is seen that the tail 
part of the neutron density originates mainly from the 
orbitals with €¢a, > —1 MeV. The average number of 
neutrons which are weakly bound or in the continuum 
is around 4.34. On the average, 2.92 of these neutrons 
are in the weakly bound orbits 3 and 4 and the others in 
the continuum. The rms radii of the core and the halo 
are 3.72 fm and 5.86 fm, respectively. A further decom- 
position shows that the two weakly bound orbitals, i.e., 
the 3rd (Q7 = 1/27) and the 4th (Q7 = 3/27), con- 
tribute mostly to the halo. This is more clearly seen in 
Fig. 3b where we represent the relative contributions of 
weakly bound and continuum orbitals to the total neu- 
tron density. The two continuum orbitals, i.e., the 6th 
(Q7 = 1/27) and the 8th (Q7 = 1/27) also contribute to 
the tail. 

If we decompose the deformed wave functions of the 
two weakly bound orbitals, i.e. the 3rd (Q7 = 1/27) 
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FIG. 3: (Color online) Neutron density distributions for **Mg. 
(a) The total density and its decomposition into core and halo. 
Contributions from several neutron orbitals around the Fermi 
level are also given. (b) Relative contributions of these neu- 
tron orbitals to the total neutron density, which is indicated 
in arbitrary units by the shaded area. 


and the 4th (Q7 = 3/27), in the spherical WS basis it 
turns out that in both cases the major part comes from 
p waves as indicated on the right hand side of Fig. 2. 
The p wave components for the 3rd and the 4th orbitals 
are 66% and 80% respectively. Having in mind that the 
occupation probabilities of these two orbitals are 0.736 
and 0.693 and each orbital has degeneracy 2, there are 
about 2 neutrons in weakly bound p states. The low 
centrifugal barrier for p waves gives rise to the formation 
of the halo. Having a small p wave component, the 6th 
orbital (Q7 = 1/27) contributes less to the halo though 
it is in the continuum and the occupation probability is 
rather large. The contribution of the 8th orbital (Q7 = 
1/27) to the tail of the density is even smaller because its 
main components are of d waves. The large centrifugal 
barrier of f states hinders strongly the spatial extension 
of the wave functions of the other two continuum orbitals, 
i.e., the 5th (Q7 = 7/27) and the 7th (Q7 = 3/27). 

In Fig. 4 the densities of the core and the halo are de- 
composed into spherical, quadrupole, and hexadecapole 
components. As is seen in Fig. 4a, the quadrupole com- 
ponent of the core is positive, thus being consistent with 


the prolate shape of 44Mg. However, for the halo, the 
quadrupole component has a negative sign, which means 
that the halo has an oblate deformation. The quadrupole 
moments of the neutron core and the halo are 160 and 
—27 fm?, respectively. This explains the decoupling be- 
tween the quadrupole deformations of the core and the 
halo as we have seen it in Figs. lb and 1c. There is 
also a noticeable hexadecapole component in the density 
distribution of the halo. 

The slightly oblate shape of the halo originates from 
the intrinsic structure of the weakly bound and con- 
tinuum orbitals. As is shown in Fig. 2, the main WS 
components of the two weakly bound orbitals, the 3rd 
(Q7 = 1/27) and the 4th (Q7 = 3/27), are p states. 
We know that the angular distribution of |Yi9(0, ¢)|? œ 
cos? 0 with a projection of the orbital angular momen- 
tum on the symmetry axis A = 0 is prolate and that 
of |Y1+1(0, ¢)|? œ sin? @ with A = 1 is oblate. It turns 
out that in the 3rd (Q7 = 1/27) orbital, both A = 0 
and A = 1 components contribute and the latter domi- 
nates. Therefore this orbital has a slightly oblate shape. 
For the 4th (Q7 = 3/27) state, there is only the A = 1 
component from the p3/2 wave, an oblate shape is also 
expected. 


In order to show that these results depend crucially on 
the single particle structure in the neighborhood of the 
Fermi surface we also investigate weakly-bound nuclei in 
the neighboring chain of Ne-isotopes. In Fig. 5a the den- 
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FIG. 4: (Color online) Decomposition of the neutron density 
of “*Mg into spherical (A = 0), quadrupole (A = 2), and 
hexadecapole (A = 4) components for (a) the core and (b) 
the halo. 
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FIG. 5: (Color online) Density distributions of ?°Ne. Details 
are given in Fig. 1 
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FIG. 6: (Color online) Density distributions of **Ne. Details 
are given in Fig. 1 


sity distributions of all protons and all neutrons in the 
prolate deformed nucleus °°Ne are shown (82 = 0.52). 
Again, as in 44Mg, due to the large neutron excess, the 
neutron density not only extends much farther in space 
but it also shows a halo structure. The neutron density 
is decomposed into the contribution of the core in Fig. 5b 
and that of the halo in Fig. 5c. In contrary to the nu- 
cleus 44Mg, we observe now a prolate halo, because the 
essential level of the halo has a large contribution from 
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the prolate A = 0 (p wave) component. In Fig. 6 we show 
similar density distributions for the oblate deformed nu- 
cleus 38Ne (82 = —0.24) which is the last nucleus within 
the neutron drip line in the present calculation. In this 
case the Fermi level is again within the pf shell. But the 
levels dominated by p waves are either less occupied or 
not so weakly bound and therefore we do not find a pro- 
nounced halo. From these examples it is clear that the 
existence and the deformation of a possible neutron halo 
depends essentially on the quantum numbers of the main 
components of the single particle orbits in the vicinity of 
the Fermi surface: s levels with A = 0 produce spherical 
halos, p levels with A = 0 prolate and p levels with A = 1 
oblate halos [8]. 

In summary, the very neutron-rich deformed nucleus 
44Mg is investigated within deformed relativistic Hartree 
Bogoliubov theory in the continuum. In contrast to sev- 
eral expectations [12, 13] a pronounced deformed neu- 
tron halo is found. It is formed by several orbitals close 
to the threshold (either weakly bound or in the contin- 
uum). They have large components of low ¢-values and 
feel therefore only a small centrifugal barrier. Although 
44Mg and its core are well deformed and prolate, the de- 
formation of the halo is slightly oblate. This implies a 
decoupling between the deformations of core and halo. 
This mechanism is investigated by studying the details 
of the neutron densities for core and halo, the single par- 
ticle levels in the canonical basis, and the decomposition 
of the halo orbitals. We also studied the weakly-bound 
nuclei in Ne isotopes and discussed the conditions for the 
occurence of a halo and its shape. It is shown that the 
existence and the deformation of a possible neutron halo 
depends essentially on the quantum numbers of the main 
components of the single particle orbits in the vicinity of 
the Fermi surface. 

Finally we note that besides the “quadrupole-bound” 
molecule [5] and the nuclear halo in deformed nuclei, 
similar coupling effects between the deformation and the 
weakly bound part of the system could also exist in other 
quantum many-body systems, such as Rydberg atoms in 
which the electron(s) can be extremely weakly bound and 
where the quadrupole moment is sizable [2]. 
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